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Abstract 

The recent identification of the dysprosium titanate compound Dy^Ti^O^ as a "Spin-Ice", i.e. the spin analog 
of regular entropic ice of Pauling, has created considerable excitement. The ability to manipulate spins using 
magnetic fields gives a unique advantage over regular ice in these systems, and has been used to study the recovery 
of entropy. Predicted magnetization plateaus have been observed, testing the underlying model consisting of a 
competition between short ranged super exchange, and long ranged dipolar interactions between spins. I discuss 
other compounds that are possibly spin ice like: Ho^Ti^Ot , and the two stannates Ho-zSn-iOi , Dy^Sn^Oi . 
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1. Introduction 

The recent realization that DTO ( -Dj^T^O?) is a 
spin analog of regular ( Ih) ice, has caused consid- 
erable excitement. Ice has fascinated several genera- 
tions of physicists in view of its apparent violation of 
the Third Law of Thermodynamics, by virtue of hav- 
ing an entropic ground state. The calorimetry experi- 
ment of Giauque and Stout [1] in 1936 was indeed one 
of the first triumphs of experimental low temperature 
physics, and theory followed experiment rapidly. The 
model of Pauling explained the origin of the entropy, as 
arising from the rearrangements of protons on the two 
possible locations on each H — O — H bond, subject 
to the Bernal Fowler ice rule of two close protons and 
two far protons for each Oxygen on the wurtzite struc- 
ture. DTO not only provides a spin realization of the 
two fold variable on each bond, it further gives one the 
handle of the Zeeman coupling between spins and the 
magnetic field, this energy reduces the degeneracy, and 
the expected entropy recovery have stimulated consid- 
erable activity in the community. In this talk I will give 
a quick review of the basic phenomena and discuss the 
modeling of the system. I will discuss the basic model 
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Fig. 1. The pyrochlore lattice 

that has been proposed by us, and tested against var- 
ious experiments, most sensitively the magnetization 
plateaus seen recently. I will also mention some unre- 
solved issues, mainly the predicted long ranged order 
in the true ground state that has evaded observation 
by neutron scattering. After summarizing the situation 
of some other candidates for spin ice behaviour, most 
notably holmium titanate and the stannates, I briefly 
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discuss the dynamical susceptibility that has been used 
to probe the nature of precursor spin liquid, i.e. the 
state that the spin ice melts into, it appears to be non- 
trivial in its correlations. 



2. Basic Phenomenon and History 

Early data [1] showed the surprising feature of en- 
tropy in the ground state, although crystallographic 
transitions prevented the cooling of Ih ice in the 
wurtzite structure to very low temperatures. Bernal 
and Fowler pointed out the importance of minimizing 
dipolar energy in the proximity of each Oxygen by 
formulating the famous "ice rule", i.e the rule that 
two protons are close by and two are further from 
each Oxygen so that there is a six fold manifold of 
states that dominate the configuration space. Pauling 
made the connection with experiments by estimating 
the global entropy arising from these local ice rules- 
the configurations on a single Oxygen influence the 
neighbours so there is a highly non trivial many body 
problem here. Pauling's estimate of the ground state 
entropy as f?ilog(|) turns out to be an inspired ap- 
proximation in that that numerical computation of the 
entropy yield estimates that are close. In the world of 
models, this led to the celebrated two dimensional ice 
model, where the combinatorial problem was solved 
exactly by Lieb[2] using Bethe's Ansatz. Anderson [3] 
showed that the identity between the Spinel B site 
sublattice and the wurtzite structure leads to an in- 
teresting connection with the configurational entropy 
in magnetite. 

Harris et al[4] showed in 1997 that the pyrochlore 
compound HTO ( Ho^Ti^O^ could be considered as 
"spin ice" by noting the absence of LRO in neutron 
scattering and from /iSR data. Moessner[5] empha- 
sized the fact that ferromagnetism plus a non collinear 
Ising type easy axis arrangement can lead to frustra- 
tion, which is commonly associated with antiferromag- 
netism. 

In a definitive experiment in 1999 Ramirez et al[6] 
showed that the entropy of DTO has a shortfall from 
the expected R log(2) by an amount that is very close to 
the Pauling value, thus establishing it as the first clear 
spin ice system. In Figure(2) we see that the entropy 
in zero field powder DTO saturates nicely, and a finite 
field, as small as .5T releases some of the entropy. We 
return to discuss the details of "entropy recovery" later. 

The currently popular basic model to describe the 
spin ice system was formulated in Siddharthan et al[7] . 
It consists of local Ising variables pointing along or 
against the local easy axis found in each case by joining 
the sites to the center of either of the two tetrahedra it 
belongs to. These local moments interact via the long 
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Fig. 2. The specific heat and integrated entropy from Rcf([6]) 

ranged magnetic dipole-dipole interaction in addition 
to the short ranged superexchange interaction J, which 
provides the single parameter of the theory. The values 
of the local moment are large for DTO (J=15/2 and 
p — gJfiB= 10 /is) so the dipolar interaction energy is 
~ 2.3°K for nearest neighbours. The dipolar interac- 
tion is long ranged but not a source of divergences in 
this problem, since all the interesting physics is away 
from zero momentum (i.e. ferromagnetic region). An 
interesting point is that the ice rule configurations are 
increasingly important at lower temperatures, in Fig- 
ure(3) we plot the fraction of tetrahedra in the ice rule 
configuration for three different values of J taken from 
a crude monte carlo sampling. It is infact easy to see 
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Fig. 3. The fraction of tetrahedra in ice ruled configurations 
for three values of superexchange J — 1.9, 1.24, .52 for curves 
from bottom up. The middle curve corresponds to DTO, and 
the other to two possible parameters for HTO. 

that the six configurations allowed by the ice rule pos- 
sess a net magnetic moment that points along or an- 
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tiparallel to the three crystallographic axes, so at low 
temperature one expects a renormalization of the ef- 
fective moment by a factor of | , and hence the uniform 
susceptibility by |. Data from Lawes and Ramirez [10] 
on susceptibility is consistent with this expectation, as 
seen in Figure(4). 
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Fig. 4. The low temperature susceptibility of HTO follows a 
Curie law but the Curie constant renormalizes to ~ ^ below 
the Ice rule scale of temperature. 

The agreement between experiments and theory on 
specific heat is quite good with the basic model[7,8], 
and one can fine-tune the basic model by adding su- 
perexchange at further neighbour distance to improve 
the agreement further [9]. We note that the most strin- 
gent test is the behaviour of the magnetization versus 
magnetic field for different directions of the field. The 
prediction of Siddharthan, Shastry and Ramirez [11] of 
magnetization plateaux for fields along < 1 1 1 > direc- 
tion represent the best test, since the plateau at ~ IT 
corresponds to some tetrahedra breaking the ice rule 
configurations, and going into a 3-in 1-out configura- 
tion. Recent experiments by Matsuhira et al[12] bear 
out this prediction very well, as seen in Figure(5). The 
onset of the second plateau is sharp at low tempera- 
tures, and can be used to experimentally define the ice 
rule energy scale. It should be regarded as the finger- 
print of the ice rule physics. This scale arises from the 
competition beween the (known) ferromagnetic dipolar 
energy and the (unknown) antiferromagnetic superex- 
change. In Figure(6) we show the range over which the 
ice rule scale changes by changing the exchange J, this 
can be used ( with linear extrapolation) to deduce J if 
we know the location of the plateau. The values of the 
magnetic moment at the first and second plateaux are 
easy to understand, the saturation values are | and | 
of the maximum. 

The problem of entropy recovery in the presence of 
a uniform magnetic field was alluded to earlier. While 
the first experiment on this in Ref[6] was on poweder 
samples, the monte carlo data [11] already gave indica- 
tions of considerable fine structure in the specific heat 
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Fig. 5. Figure on top is the theoretical prcdiction[ll], the 
dotted curve should be compared with the one below from 
experiment [12] . 



T=0.4K, [111] direction 




- J_s = 0.52 
J_.s=l 24 
I_s=1.92 



Fig. 6. Three values of J and the resulting range of the plateau, 
the values of J are 1.92, 1.24, .52 from left to right. 

depending upon the direction of the field relative to the 
symmetry axes. For example a double peak structure 
was predicted for fields along < 1,-1, V2 > direction 
which looks rather close to the double peak seen very 
recently in experiments[12]. 

We next mention holmium titanate HTO, where 
Ho 3+ has an identical magnetic moment as Dy +3 , 
and hence the same dipolar interaction, but possibly 
different exchange J. Initial experiments [6] showed a 
rising specific heat down at ~ .5°K, at which point the 
system fell out of equilibrium. One can interpret this 
in one of two possible ways. It is imaginable [7,9,11] 
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that the system undergoes a transition to an long 
range ordered (LRO) state which is anyway predicted 
by theory (see below), and that the slow dynamics 
hides this transition. One has to then explain the dif- 
ference between DTO and HTO as possibly arising 
from differences in the ice rule scale, or equivalently 
J, so that HTO would have a much lower ice rule 
temperature ( hence larger |J|) whereas DTO with a 
higher ice rule temperature would be stuck in a subset 
of configurations that would prevent it from under- 
going a transition to a state with lower free energy. 
The other view point [14] ascribes the difference in 
behaviour to the significant hyperfine coupling con- 
stants in HTO, as known from early work of Blote et 
al[13]. Subtracting the nuclear component makes the 
data for HTO look similar to DTO, and the recent 
first principles calculations of the hyperfine constants 
of these two compounds [16] lends weight to this view 
point, as does the absence of LRO as indicated by 
neutron scattering [14,15]and other experiments at 
low temperatures ~ 50° mK. 

Staying with the issue of LRO, we mention that even 
in the case of DTO one expects LRO at low enough 
temeratures. Siddharthan et al [9,11] were the first to 
show, by explicit enumeration of ground states for finite 
clusters, that the model for DTO should have LRO of 
a certain type. This has been corraborated by Melko 
et al[17], who use a loop algorithm to equilibriate the 
system rather than single spin flips, and find the same 
structure, with a transition at ~ .2°K. In contrast, 
neutron scattering[18] sees no signs of LRO down to 
low temperatures ~ 50°rnK ! Therefore it seems to me 
that DTO and possibly HTO may yet have a surprise 
in store at low temperatures. 

Turning to other possible spin ice compounds, the 
stannates M2S712O7 with M = Dy,Ho have been re- 
cently proposed [19]. These have small differences in 
lattice constants, and hence slightly different dipolar 
as well as superexchange interactions, and should be 
useful in fixing several details of the models. 

Finally on the topic of spin ice, I would like to men- 
tion two very recent and nice experiments on the AC 
susceptibility of DTO. Schiffer et al [20] and Matsuhira 
et al [21] have found that the curie type divergence of 
the low field susceptibility ( as in Figure (4)) is infact 
cut off at a low temperature ~ 15° K that depends 
upon the frequency of the probe field, resulting in two 
maxima, one at a higher temperature ~ 18° K and the 
other at a lower temperature ~ 1.5° K. Each of these 
maxima is freqency dependent, and the from the Ar- 
rhenius dependence of these we can extract a character- 
istic energy scale e c corresponding to these peaks. The 
higher one yields an energy scale ~ 220° K, very nearly 
the crystal field splitting seen in [7], and the lower one 
give ~ 10° A", presumably a multiple of the ice rule en- 
ergy scale. The real as well as imaginary parts of the 



susceptibility show the same features, and these seem 
to be consistent above ~ 10° K with a dynamical sus- 
ceptibility n^torT" ' w here \o ~ y and a departing 
from the Debye value of unity, and further being rather 
temperature dependent, with a ~ .544 at 17°K and 
da/dT ~ .06/° K. Thus the non Debye relaxation with 
temperature dependent a, and the origin of the relax- 
ation rate r ~ to exp(e c /fcT) demand a fundamental 
understanding that is missing at the moment. 
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